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Abstract There are various types of global and local spacetime invariant in general rel¬ 
ativity. Here I focus on the local invariants obtainable from the curvature tensor and its 
derivatives. The number of such invariants at each order of differentiation that are alge¬ 
braically independent will be discussed. There is no universally valid choice of a minimal 
set. The number in a complete set will also be discussed. The invariants can then be used to 
characterize solutions of the Einstein equations (locally), to test apparently distinct solutions 
for equivalence, and to construct solutions. Other applications concern limits of families of 
spacetimes, and the characterization of horizons and singularities. Further uses are briefly 
mentioned. 
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1 Introduction and motivation 

This lecture developed from one given some years ago at the retirement party of my friend, 
colleague and co-author Dr. Eduard (“Eddie”) Herlt. At the Lahore meeting, time did not 
permit a complete exposition of all aspects of the subject. This text follows what was said 
rather than what might have been said in an expanded version. I intend to write up the 
expanded version as a review article in due course. 

I could have been more precise in my title, at the cost of being rather lengthy. There are 
many occurrences of the word ’invariant’ in our field, e.g.: 
gauge-invariant (in gauge theory) 
gauge-invariant (in perturbation theory) 
invariant under a symmetry 
Lorentz-invariant 
scale-invariant. 

And I could be talking about some global conserved quantity, without a well-defined local 
density, e.g. Bondi mass. 
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What I actually meant was local, geometric, invariants of spacetime: essentially curva¬ 
ture invariants. 

I started to be seriously interested in this area when we were writing the first edition 
of the exact solutions book (II] in the late 1970s. One of the big problems we faced was 
that of identifying the same solution when found with different assumptions or for different 
reasons, and presented in different coordinate systems. I realised at the time that the work 
of the Stockholm group on invariant classification, which I will mention later, and which I 
first heard about in the late 70s, held the key. 

This first application, as far as I was concerned, the “equivalence problem”, was what 
got me started. 


2 Mathematics 

2.1 Definitions of invariants 

Christoffel proved in 1869 that scalars constructed from the metric and its derivatives must 
be functions of the metric itself and the Riemann tensor and its covariant derivatives. 

The first examples to spring to mind are scalar polynomial (s.p.) invariants, such as 
or CabcdC‘^‘^"^Cef‘‘^■ Often when people just say “invariant” they mean “s.p. invari¬ 
ant”. 

However, these are not adequate in all circumstances, as one can see by noting that pp 
waves and flat space both have all scalar polynomial invariants, of all orders, equal to zero 

i). 

Fortunately they are not the only choice. An important alternative is to use ideas due to 
Cartan, as follows. As a side benefit, in my view an important one, Cartan invariants require 
less calculation, in general. 

Let F(^) denote a “suitable” frame bundle over a spacetime ^ (i.e. take the set of all 
frames at each point) and be the set {Rabcd , Rabcd j , • • ■, Rabcd Ji A ■■■/,,} of the components 
of the Riemann tensor and its derivatives up to the qth in a frame. 

Choose from in a canonical and invariant way, e.g. use the principal null direc¬ 

tions of the Weyl tensor, when they are distinct, as the basis vectors. The resulting 

{Rabcd ^ Rabcd'.f ^ • ■ ■ ^Rabcd;f[f2---fq{ 

are called the Cartan invariants. They are scalars, because the frames are invariantly de¬ 
fined, e.g. if a, b, c and d are basis vectors (not necessarily all distinct) of the chosen frame, 
one of the Cartan invariants is Ripidy. The idea is like characterizing a symmetric 
bilinear map (matrix) by its eigenvalues. 


2.2 Independence of (s.p.) invariants 

In a manifold of n dimensions, at most n scalar invariants can be functionally indepen¬ 
dent, i.e. independent functions on ^. 

The number of algebraically independent scalar polynomial invariants, i.e. s.p. invari¬ 
ants not satisfying any polynomial relation (called a syzygy), is rather larger, as we shall see 
next. (A set of algebraically independent Cartan invariants in a general spacetime, written 
using the Newman-Penrose spinor formalism, is given in it takes fully into account the 
Ricci and Bianchi identities.) 
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Larger still is the size of a complete set of s.p. invariants (a finite Hilbert basis): such a 
set {/i, / 2 ,..., 4} is complete if any other such s.p. invariant can be written as a polynomial 
in the Ij but no invariant in the set can be so expressed in terms of the others. 

One way to find the number of algebraically independent scalar invariants is to consider 
Taylor expansions of the metric and of the possible coordinate transformations 0]; another 
way follows Hilbert 0. 

The result 0 is that in a general V„ the number of algebraically independent scalars 
constructible from the metric and its derivatives up to the pth order (the Riemann tensor and 
its derivatives up to the {p — 2)th) is 0 for p = 0 or p = 1 and 


N{n,p) 


«[«+ l][(n + p)!] 

2«! p\ 


(w + p+1)! 

(n-l)!(p+l)!^”’ 


for p > 2, except for A^(2,2) = 1. Thus in a general space-time the Riemann tensor has 
A^(4,2) = 14 algebraically independent scalar invariants. In particular cases the number is 
reduced. 

The origin of many syzygies can be understood in terms of the vanishing of anyobiect 
skewed over (n-|- 1) indices in n dimensions SSi- In a series of papers iT(iTTll fI3l 
Carminati and Lim have given a detailed discussion of the syzygies for scalar polynomial 
(s.p.) invariants of the Riemann tensor, using graph-theoretic techniques. 

A given invariant may be written in more than one way, due to symmetries, and other 
relations between components, of the Riemann tensor and its derivatives. This is essentially 
a problem in representations of the permutation group. The issue is to select a canonical 
representative of each equivalence class in the orbit under permutations. 

One wants to select a canonical set of invariants and then express any other invariant in 
terms of canonized members of the algebraically independent set. 

Several methods have been used to do this for s.p. invariants, e.g. by Homfeldt, by 
Fulling et ah, by Ilyin and Kryukov, and by Dresse. The method most readily available is 
due to Portugal It has been implemented in Maple™ and Mathematica® by Martin- 
Garcia et al, and is distributed in xAct, a package for use with Mathernatica®. It has (e.g.) 
been applied to all objects with up to 12 derivatives of the metric [3l. 

Any complete set of s.p. invariants of the Riemann tensor, and any set which always 
contains a maximal set of independent scalars, contains redundant elements. Hence all the 
old papers giving specihc sets of 14, in 4 dimensions, are inadequate. 

The smallest known set of s.p. invariants for a general Riemann tensor that always con¬ 
tains a maximal set of algebraically independent scalars consists of 17 polynomials Q, 
though 16 suffice for perfect fluids and Einstein-Maxwell fields. There is a special subclass 
of spacetimes (ii which require 18. 

The smallest set of s.p. invariants known to be complete contains 38 scalars IT^ . and 
Lim and Carminati in proved that it is minimal. 


2.3 Computation of invariants 

The s.p. invariants are expressions whose contractions hide very large numbers of individual 
terms, and are therefore hard to calculate. It is thus very useful to adopt a method that reduces 
the number of terms. One can use a bivector formalism (see e.g. IT3l l. Another way, in 
4 dimensions, is to use the Newman-Penrose complex spinor (or null tetrad) formalism, 
the GHP formalism, or related ideas: this has analogues in other dimensions (e.g. in three 
dimensions one can use real two-component spinors: see d). 
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For higher dimensions it may be efficient to use parallel computing, i.e. to split the 
calculation into subproblems e.g. by fixing some of the dummy indices, and then sending 
the subproblems to separate processors 11911 . 

A simple practical way in many cases is to use Cartan invariants instead. They are quite 
cheap to calculate provided the frame choice is easy to calculate with, in particular because 
one never has to multiply curvature components, or their derivatives, together. In a general 
case, either s.p. or Cartan invariants would characterize the spacetime, as we shall see below, 
so the choice is really a question of convenience or purpose. There is an open research 
problem about which is actually more computationally efficient (or for which cases). 

In four dimensions we can always choose a frame, for example by the principal null 
directions of the Weyl tensor, and thus compute Cartan invariants, but in general this may 
entail using the unpleasant algebraic expressions which arise from the general formulae for 
solutions of quartics. However, in practice many solutions give way easily. Note that in 5 
dimensions or more there is no guarantee we can calculate frames of a specified type even 
if we can show they must exist, as there is no general algebraic formula for solutions of 
quintics. 


2.4 Do scalar polynomial invariants suffice? 

Until 2009, 1 would have said definitely not. I mentioned earlier that pp waves and flat space 
both have all scalar polynomial invariants, of all orders, equal to zero. In fact all vacuum type 
N and III metrics with p = 0 have this property lioll . There are also metrics which have equal 
non-zero s.p. invariants e.g. l2lLl22Ll23ll . 

These ambiguities are associated with the indefiniteness of the metric and the non¬ 
compactness of the Lorentz group f^ . 

Coley et al. gave an argument that spacetimes are completely characterized by their 
s.p. invariants, except for spacetimes in the Kundt class. Although their proof seems to me to 
have a gap, and I have not yet been able to see if it has been filled, the result seems correct. 

In particular Coley and collaborators have followed up by a substantial number of pa¬ 
pers, which I did not have time to review comprehensively in my Lahore talk. These consider 
not only 4D spacetimes, but higher dimensional spaces and various signatures. I just men¬ 
tioned a few to give the flavour. In the initial paper 12^ . they also discuss what properties a 
given set of invariants characterize. 

- Hervik and Coley show that any metric with an analytic continuation to the pure 

Riemannian case is completely characterized by its s.p. invariants 

- Coley et al. give a detailed analysis of the generic 5D spacetimes. 

- Hervik gave arguments that spacetimes not characterized by s.p. invariants are lim¬ 
its of families that are so characterized. 

It may be worth noting that since it is obvious that the algebraically independent set of 
Cartan invariants given in completely and uniquely specify the s.p. invariants, the result 
of i^l amounts to saying that in general the converse is true, i.e. that the s.p. invariants 
uniquely determine the Cartan invariants (given the procedure for the choice of frame). 

The work just described picked out the Kundt class as the exceptions. These are space- 
times which have a null vector held which is geodesic, and expansion-, shear- and twist- 
free. In a paper which appeared just before this class was discussed in detail. For 
it, one has to use the Cartan approach or some equivalent. 
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In 12^ Coley et al. discussed the Petrov D examples in the Kundt class. Podolsky and 
Svarc 13011 gave a classification of Kundt metrics in arbitrary dimension, and later discussed 
their physical interpretation. A subclass of “universal spacetimes” was discussed by Hervik 
et al. 13111 : universal in that they solve vacuum equations of all gravitational theories with 
Lagrangian constructed from the metric, the Riemann tensor and its derivatives of arbitrary 
order. Hervik et al. found all cases where there exists a continuous family of metrics 
having identical [scalar] polynomial invariants. 


3 Uses of spacetime invariants 

3.1 The Equivalence Problem 

Now let me return to my original motivation, i.e. identifying the same solution when found 
with different motivations and assumptions and presented in different coordinate systems. It 
could be considered either a mathematical or a physical problem. It is essentially resolved 
by the following theorem (in terms of Cartan invariants). 

Theorem 1 (Sternberg 1^ . Ehlers 1^ ) Given two spacetimes, ^ and each expressed 
using some frames E and E, then there is an isometry which maps (x, E) to {x, E) if and 
only if for .M is such that the independent quantities and the functional relations are 
the same for and , where p is the last derivative at which a new functionally 

independent entry in arises. 

See for a more careful statement. 

One could therefore calculate everything on a suitable frame bundle. 

- Christoffel used coordinate frames. This works but the dimension of F is large. 

- Cartan proposed using frames with constant scalar products, e.g. null frames or or¬ 
thonormal frames, which is better (hence the credit for Cartan invariants). 

- In 1965 Brans proposed practical use of this method via lexicographic indexing of 
components. 

- Later Brans and, more fully, Karlhede realised that a more efficient way to im¬ 

plement the idea was to restrict the frames wherever possible by only allowing canonical 
forms of the curvature and its derivatives (i.e. go to Cartan invariants). 

- This makes the whole thing manageable (with a computer) and it was first implemented 
by Aman in 1979 

Given the Coley et al. results, we now know we could in general work only with s.p. 
invariants, but the present software SEl] uses the Cartan method. In an algebraically 
general spacetime, one can choose the principal null directions of the Weyl tensor to fix the 
frame, or in a conformally flat spacetime one can begin with eigenvectors of the Ricci tensor. 
In the actual implementaion we use the Weyl PNDs, where they exist, as the first choice. 

Many details have been considered (by Aman, myself, Skea, d’lnverno, McLenaghan, 
Pollney and others) in specifying frames and in making the computations more efficient. 
The main aspects are the enumeration of canonical forms, the tests of which canonical form 
applies, and the transformation of the non-canonical to the canonical. The literature cited in 
the exact solutions book gives more information. 

A first extension is to maps more general than isometries, e.g. homotheties 14^] or con¬ 
formal equivalence 1^ . One can of course tackle Euclidean metrics and metrics in any 
dimension, though the details of making the method work efficiently are different (e.g. I^ I. 
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The above ideas are in fact a special case of Cartan’s general procedures (see e.g. (4^ 1. 
which apply to other situations which can be expressed in differential forms and a connection 
or similar structures. In particular it applies to the equivalence of (systems of) differential 
equations, under coordinate transformations (see the books of Gardner and Olver El). 

Another context is that of gauge theories of physics, the gauge potential being just a 
connection lE^ . 

One interesting issue for spacetimes is how large p in the theorem has to be. Karlhede 
ifssll showed p <1. For a long time examples suggested p < 3. At the time of the second 
edition of the exact solutions book an example with p = 5 was known, as were a number 
of bounds for subclasses. We quoted results suggesting that p <5 was generally true. 

In 2009 (the result was actually hrst announced in 2007), Milsom and Pelavas 
showed that the original bound of 7 is sharp, by giving a (unique) example. This is still 
small enough for calculation to be practical. The example is of course in Kundt’s class, as 
are those further examples they found where p = 6. 

As well as giving a way to compare two spacetimes, the theorem above also implies 
that the Cartan invariants uniquely classify and determine the local geometry. This implies 
that all local properties are encoded in this information. I will rather briefly discuss some 
of the possible resulting applications. Aspects I was not able to cover in Lahore included: 
applications to three-dimensional spaces, incuding the use in checking junction conditions; 
applications in numerical relativity; applications to cosmology and to perturbations; use 
in constructing alternative theories of gravity; and applications to internal state spaces of 
thermodynamics. I think there must be more we have not yet explored. 

Note that global topology is of course not determined by the invariants we have been 
considering (local invariants cannot distinguish a plane and a flat torus) and nor are contin¬ 
uations which do not respect the required differentiability for the theorem above. 

My original motivation can now be carried out for many cases. I used these methods, for 
example, to disentangle the known exact solutions for inhomogeneous cosmologies with a 
symmetry group G 2 . It would be good to have a complete online catalogue of (at least) the 
solutions in 13511 . Skea and Lake have independently implemented systems to facilitate this. 
Properties such as the isometry group can also be found l50l l5lL 15211 : for examples see e.g. 

Similar results apply for homotheties E^l. 


3.2 Finding solutions 

The same ideas can be used to find solutions, by working out consistent sets of Cartan 
invariants and then integrating to get a metric. Examples are mentioned in the exact solutions 
book. Machado Ramos and Edgar used these ideas, implemented in their invariant 

operator formalism, to hnd pure radiation solutions of types N and O with a cosmological 
constant. 

Coley and collaborators looked for all spacetimes in which all s.p. invariants vanish 
(which they called VSI spacetimes although not all Cartan scalar invariants vanish). In the 
4-d Lorentzian case these give all spacetimes indistinguishable from flat space by s.p. in¬ 
variants i^ . The higher dimensional cases give examples in supergravity and string theory. 
A closely related series of papers has studied the spacetimes with constant s.p. invariants 
(CSI spacetimes). For a review of this work see 1^ . 
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3.3 Limits of families of spacetimes 

Another aspect, described first by Paiva et al. is the use of invariants to work out the 
limits of families described by parameters. This enables one to find all possible limits in 
a coordinate-free way, whereas previous treatments were in effect trial and error. This is 
another area where further work would be useful. 

I did not give more detail in my Lahore lecture, but passed quickly on to limiting points 
within spacetimes. 


3.4 Black hole and other horizons 

Karlhede et al Iboll first noted that Rabcd\eR‘‘^‘^'^'‘ = 0 at the Schwarzschild horizon (so a 
prudent space traveller might monitor that). Skea in his thesis (HI] noted that this is not true 
for other horizons (a point rediscussed by Saa for higher-dimensional static cases: Saa 
also found points where Rabcd\eR“^‘^‘^'‘ = 0 but which are not horizons). Lake l[^ contin¬ 
ued the work on Kerr by considering hrst derivative invariants, and found their vanishing 
characterized the horizons. 

Related to my arguments against the claims of Antoci and others that the Schwarzschild 
horizon was singular, I proposed in 2006 a new test for occurrence of a horizon using 
ratios of Cartan invariants which works in all cases of Petrov type D. I believe it works more 
generally, but this needs further study. Maybe it could also be useful to numerical relativists 
(and space travellers). 

More recently Abdelqader and Lake 1^ . have given an invariant characterization 
of the Kerr horizon, and, prompted by that. Page and Shoom have given another for 
stationary black holes. I have not yet had the opportunity to check how these relate to earlier 
characterizations. 

Moffat and Toth considered the relation of the “Karlhede invariant” (i.e. Rabcd-,eR‘'^‘^‘^’'^ 
to discussions of a “hrewall” at the horizon. 


3.5 Singularities 

Singularities in general relativity are defined to occur when a causal geodesic cannot be 
continued to inhnite affine parameter values even when the spacetime is maximally extended 
- “geodesic incompleteness” (for the reasons for this dehnition see (b^). 

It is well-known that this happens if an s.p. invariant of the Riemann tensor itself (not 
its derivatives) blows up along the geodesic. This does not however mean that: 

(i) an “infinite” Riemann tensor implies a singularity, or 

(ii) if an invariant blows up, there is a singularity, or 

(iii) at all singularities, an s.p. invariant of the Riemann tensor blows up. 

1. An “infinite” Riemann tensor does not imply a singularity. Geodesics can be continued 
across a delta function curvature modelling a thin shell or an impulsive gravitational 
wave. 

2. The blowing up of an invariant does not imply a singularity. For example the in¬ 
variant I /Rabcd-,eR“^‘^'^’‘ blows up at the Schwarzschild horizon, but the horizon is not 
singular. 
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3. There are singularities at which no s.p. invariant of the Riemann tensor blows up. At 
“whimper” singularities an invariant involving first derivatives does, though An 
example was studied by Podolsky and Belan 17IH . 

Question: when does blow up of higher derivative invariants imply a singularity? 
Another application is to “directional” singularities, where a singular point apparently 
has directionally dependent limits. Scott and Szekeres f^l?^ showed that the directional 
singularity of the Curzon metric hid more extended regions at whose boundary the original 
coordinates broke down. My student Taylor i74ll showed that such cases could be appropri¬ 
ately “unravelled” by using level surfaces of Cartan invariants to define new coordinates. 
Lake used the Weyl tensor s.p. invariants to show that the Kerr singularity was not 
directional. 

With Coley and others iIt^ . I considered “kinematic singularities” (where fluid expan¬ 
sion blows up), giving examples in which, given an integer p, the Cartan (and hence s.p.) 
scalars can be finite up to the p-th derivative, but not the {p + l)-th. 

Geodesic continuation needs a metric. In invariantly-defined frames the connection 
coefficients (which would be C*^) are typically expressible as ratios of first derivative Cartan 
invariants to zeroth derivative ones. We know that there are “intermediate” or “whimper” 
singularities where s.p. invariants of the Riemann tensor do not blow up, while s.p. invariants 
of the first derivatives of the Riemann tensor do. Hence I conjecture that under some suitable 
differentiability conditions: 

Conjecture: Spacetime singularities are either locally extendible or at least one Cartan 
invariant in has an infinite limit along any curve approaching the singularity. 
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